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We perform a non-perturbative determination of the 0(a)-improvement coefficient csw for the 
^^ ' Wilson quark action in three-flavor QCD with the plaquette gauge action. Numerical simulations 

^ I are carried out in a range of /3 = 12.0 - 5.2 on a single lattice size of 8'^ x 16 employing the Schrodinger 

OO ■ functional setup of lattice QCD. As our main result, we obtain an interpolation formula for csw and 

Cn ' the critical hopping parameter Kc as a function of the bare coupling. This enables us to remove 

^^ ' 0{a) scaling violation from physical observables in future numerical simulation in the wide range of 

>0 , /3. Our analysis with a perturbatively modifled improvement condition for csw suggests that flnite 

^—f^ ' volume effects in csw are not large on the 8^ x 16 lattice. We investigate Nf dependence of csw 

by additional simulations for Nf = 4, 2 and at /3 = 9.6. As a preparatory step for this study, we 
also determine csw in two-ffavor QCD at /3 = 5.2. At this /3, several groups carried out large-scale 
calculations of the hadron spectrum, while no systematic determination of csw has been performed. 



^ ■ I. INTRODUCTION 



Recent large-scale simulations in two-flavor QCD have demonstrated that quenching artifacts found in various 
, . , physical observables are reduced by dynamical effects of up and down quarks. There has been significant progress 
j^ ' also in the algorithms for QCD with odd numbers of flavors: while the conventional Hybrid Monte Carlo (HMC) 
algorithm flj is applicable only to even-flavor QCD, the exact algorithms, such as the Multi-Boson ^ and the 
Polynomial HMC algorithms yf capable of odd-flavor cases, have been developed. Clearly the time has come 
to undertake a fully realistic and extensive simulations of QCD with dynamical up, down and strange quarks. 

Since simulations with dynamical quarks are computationally demanding, highly improved lattice actions 
should be employed in the three-flavor simulations. The leading cutoff effect in physical quantities is 0{a) with 
the Wilson quark action, and this error can be removed by adding a single counter term, the Shcikholcslami- 
Wohlert (SW) term Q, to the action with non-perturbatively determined coefficient csw- However, csw has 
been determined only in quenched and two-flavor QCD so far 



In this article, we perform a non-perturbative determination of csw in three-flavor QCD with the plaquette 
gauge and the Wilson quark actions. In Rcfs. [3,1a]' however, we found that this theory has a non-trivial phase 
structure: there is an unphysical phase transition at /3«5.0, where the lattice cutoff a~^ is roughly 2.6 GeV. It 
is expected that the phase transition strongly distorts scaling properties of physical observables. This transition 
is considered as an artifact due to the finite lattice spacing and can be removed by the use of improved gauge 
actions [3,|3]- Therefore, there are two strategies for meaningful simulations in three-flavor QCD: i) use the 
plaquette gauge action at /? ;» 5.0, or ii) use an improved gauge action, if simulations at a^^ < 2.6 GeV 
are needed. We explore the former strategy in this article, and the latter possibility is studied in a separate 
publication |^. 

In our determination of csw> we follow the method proposed by the ALPHA Collaboration in Refs. Ia,ll3l- 
We explore a wide range of /3= 12.0-5.2, which is significantly higher than the phase transition point /3«5.0, 
employing a single lattice size of 8^ x 16. As our main result, we derive an interpolation formula for csw as a 
function of the bare coupling. The critical hopping parameter Kc in the 0(a)-iniproved theory is also obtained 



as a byproduct. We examine finite volume effects in csw by modifying the improvement condition at one-loop 
order of perturbation theory. Additional simulations in four-, two-flavor and quenched QCD at /3 = 9.6 are 
carried out to investigate the Nf dependence of csw- 

As a preparatory step toward this study, we also determine cgvv in two-flavor QCD at f3 — 5.2. In the 
previous work by the ALPHA Collaboration 6], they carried out the non-perturbative tuning of csw at /3>5.4, 
and derived an interpolation formula of their result as a function of ^q. However, due to the limitation of 
available computer power, recent large-scale simulations by the UKQCD 11] and JLQCD Collaborations [l^ 
were performed at a lower value /3 — 5.2 with csw obtained by extrapolating the ALPHA'S formula. We, 
therefore, determine cgw directly at P — 5.2 in order to see if the extrapolation of the formula to this value of (3 
really works out, and to confirm that 0{a) scaling violation is absent in the UKQCD and JLQCD simulations. 

This paper is organized as follows. In Sec^ we briefly introduce the method for the non-pertubative tuning 
of Cgw employed in this study. Section lTTTl is devoted to detailed description of our numerical analysis and results 
in two-flavor QCD at /3 = 5.2. We present our results in three-flavor QCD and discuss its 0{a/L) uncertainty 
in Sec. IIVI Finally our conclusion is given in Sec. 



II. IMPROVEMENT CONDITION FOR 0(a)-IMPROVEMENT 

In our determination of cswj we basically follow the method proposed in Refs. |a. Il0l|. which employs the 
Schrodingcr functional (SF) setup of lattice QCD [l^. In this section, we briefly introduce the SF setup and 
the choice of the improvement condition to flx csw- 



A. SF setup 

The SF is the generating functional of the fleld theory with the Dirichlet boundary condition imposed in the 
temporal direction. In this study, the spatial link variables at the boundaries are set to the following diagonal, 
constant SU{3) matrices 

Uk{x,XQ)\^^^exp[aCk], C/^ (x, a;o)|^g^j. = exp [aC^] , (1) 

(2) 





where Lk {k — 1, 2, 3) and T are physical lattice sizes in the spatial and temporal directions. All quark fields at 
the boundaries are set to zero. In the spatial directions, the periodic boundary condition is imposed for both 
gauge and quark fields. 



We use the plaquette gauge action 



^9 = f E Tr [1 - ^../.^] , (3) 
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where Ux.fj,u is the product of gauge link variables Ux,fi around the plaquette 
The 0(a)-improved Wilson quark action ^j is given by 

Sq = ^QxDxyqy, (5) 



x,y 



Dxy = 5xy - KyUl- 7^) Ux,fiSx+f:i,y + (1 +7m) Ul-u.u^x~fi-:v\ + ■^KcswQ-tJ.uFx.p.^Sxy, (6) 



with the field strength tensor Fx^^i, defined by 



^x^fiu — g xy^x.^ii' I" ^x^u—p- ' ^ x, — f.i — u ^" ^x, — u^) y^'^'ji t v' / 



where (h.c.) denotes the hermitian conjugate of the preceding bracket, and cr^i/ = (i/2) [7^,7,^]. The last term 
in Eq. © is the counter term to remove 0{a) effects in on-shell quantities. Its coefficient csw is set to unity 
to remove tree-level 0{a) scaling violation from physical observables. The main purpose in this article is non- 
perturbative tuning of csw for removal of all 0(agQ ) scaling violation (n > 0). For the 0(a)-iinprovement of 
the SF itself, we add counter terms made of the gauge and quark fields at boundaries to the lattice action. 
However, these counter terms affect the PCAC relation at order of a^ or higher, and hence are not necessary 
for determination of csw from the PCAC relation. In the ALPHA Collaboration's studies, the counter terms 
are omitted except for a term 

SSg = ^ (c, - 1) Y. Tr [1 - C/,,^o] , (8) 

x,/i {xQ—0,T—a) 

which is made of the temporal plaquettes touching the boundaries. In this study, we also include this counter 
term to the total lattice action S = Sg+dSg+Sq so that we can directly compare our and ALPHA Collaboration's 
results. The coefficient of the counter term ct is set to the one- loop estimate in Ref. J14l |. 



B. Improvement condition 

We determine csw by imposing the validity of the PCAC relation 

i(9^ + a;)AC„p,^ = 2mF^ (9) 

up to order of a^. The pseudo-scalar operator and 0(a)-improved and unimproved axial currents are given by 

P" = V^Tsr^, (10) 

^Snp,p - ^;; + CA^(a^ + 9;)p", (11) 

A''^ = ^7M75r>, (12) 

where d^ and d* are the forward and backward lattice derivatives and SU{Nf) generators r" act on the flavor 
indices of the quark fields 'tp and -0. 

We measure two correlation functions 

fAixo) = -^^^(AS(x)O'^), (13) 

fpM = -^^^(P"(x)0"), (14) 

where (• • • ) denotes the expectation value after taking trace over color and spinor indices and summing over 
spatial coordinate x. For the source operator, we take 

O" = a6 5]C(y)75r"C(z), (15) 



y>z 



defined from the boundary fields 



^W = ^> ^"« = 4^' ^''^ 



where p(x) is the quark field at xg = and is set to zero in the calculation of /^ and fp. The bare quark mass 
is then calculated from /a and fp through the PCAC relation Eq. ^■. 

■m{xo) = r{xo) + CAs{xo) (17) 

r{xo) = \{do + d*o) fAixo)/ fp{xo) (18) 

s{xo) = -^adQdQfp{xo)/fp{xo). (19) 



We can calculate another set of to', r' and s' from the correlation functions 
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f'AT-xo) = +^^:^^{A^^{x)0'n, (20) 



Nj-l 
fpiT-xo) = -jJ-^{p-{x)0'n, (21) 

using the source operator at the other boundary 

O'''^ = a6^C'(y)75r'^C'(z), (22) 

where C is the boundary field at xq = T. 

The improvement condition to fix csw is obtained by requiring that quark masses calculated with different 
boundary conditions coincide with each other. However, a naive condition to(xo) = m'{xo) requires a non- 
perturbative tuning of ca as well as csw- To eliminate ca from the process, it was proposed in Ref. 5' to use 
a modified definition of the quark mass 

M(xo,yo) = m[xo) — -—^s{xo), (23) 

and similarly defined M'{xo,yo)- Therefore, csvv is tuned so that the following mass difference 

AM{xo, 2/0) = Mixo, yo) - M'{xo, yo) (24) 

vanishes with a certain choice of {xo, yo)- 

In principle, we can take an arbitrary choice for {xo, yo), since a change of the choice leads to a difference in 
0(0^) scaling violation in physical observables. In this study, we take {xo,yo) — (3r/4,T/4) for AM(a;o,?;o), 
and {xo,yo) — {T/2,T/4:) for M(a:o,yo)- The latter is used to specify the massless point. We note that this 
choice is the same as that in the ALPHA'S studies in quenched and two-flavor QCD. From now on, M and AM 
without arguments denote M{T/2,T/i) and AA/(3r/4,T/4), respectively. 

In practice, cgw is determined by demanding that AI and AM satisfy the following improvement condition 

M = 0, 



AM = AM(o), 



(25) 



where AM'^*'^ is the tree-level value of AM at the massless point M = on the finite lattice volume L"^ x T. We 
tune AM to AM^*^' but not to zero so that the weak coupling limit of the non-perturbatively determined csvv 
is exactly unity. On our lattice size of 8^ x 16, aAAf(°) =0.000277 5J. We also note that the tuning of M to the 
massless point provides a non-perturbative estimate of the critical hopping parameter Kc in the 0(a)-improved 
theory. 



III. TWO-FLAVOR QCD AT /3 = 5.2 
A. Simulation method 

In this section, we report the determination of cgw in two-flavor QCD at /3 = 5.2. 

Our numerical simulations are carried out on a S'^ x 16 lattice at six values of csw in a. range csw = 1-5 -3.0. 
We choose two to four values for the hopping parameter K at each csw so that we have data of AM at both 
of positive and negative values of M, and/or at M close to the massless point {\aM\ <0.01 in our study). This 
enables us to tune (M, AM) to (0, AM*^"^) by an interpolation or short extrapolation. The simulated values of 
Csw and K are summarized in Tabled 

We use the standard HMC algorithm with the asymmetric even-odd preconditioning described in Refs. [isl 
[la, [13 foi' the determinant of the quark matrix D. We solve the linear equation DX = B using the BiCGStab 
algorithm '18'| with the stopping condition 

1-14 



Rm\B\\<lO-'\ (26) 



where Ri = DXi — B is the residual vector and Xi is the estimate for the solution X in the z-th BiCGStab 
iteration. The HMC trajectory length is fixed to the unit length. We set the number of the Molecular Dynamics 
steps to 60-80, which achieves the acceptance rate higher than 80%. 

After the thermalization of 500 HMC trajectories, we accumulate the statistics A^traj summarized in Table Q] 
The correlators fx and f'x {X = A, P) are measured at every trajectory. We use the jackknife method to 
estimate statistical errors of fx , f'x and all results derived from them. 



B. numerical result 

Numerical results of ill and AM are summarized in Table ^] In order to fix csw and Kc satisfying the 
improvement condition Eq. H25|) . we parameterize M and AM by a simultaneous fit in terms of 1/ K and csw'- 

o\J b\„' (^\ (n\ n dn 



aM = aM+"-^ + -^ + c^ csw + 4'^ ^sw' + ^ cgw, (27) 

(1) f,(2) 



K K^ ^* ^'^ K 



aAM = aAM + %^ + ^ + 4m ^sw + 4m ^sw' + ^ csw (28) 

Fit parameters are summarized in Tablc lTlTl Figure^shows 1/K dependence of M and AM, and AI dependence 
of AM at each csw- By interpolating (A/, AAf ) to (0, AAi^'^^) with this parameterization, we obtain 

csw = 1.908(64), (29) 

Kc = 0.1381(12). (30) 

We also test another method for the parametrization of M and AM in order to estimate the systematic error 
due to the simultaneous fit Eqs. H27|l and (|28(l . At each cswi we determine AM at M = by a linear fit 

aAM = a'„i + b'^iaM). (31) 

At Csw = 1-90 and 2.02 where we simulate more than two values for K, we also test a quadratic form and find 
that the higher order contribution is small and can be safely neglected in this analysis. Figure El shows csw 
dependence of AM at M = which we parametrize by linear or quadratic forms 

aAM = a^J^J + b'^^^ csw + c^^/ csw^- (32) 

By tuning AM to its tree-level value, we obtain csw — 1.979(68) from the linear fit as reported in Ref. (ia |. 
and 1.975(50) from the quadratic one. These are consistent with the result from the combined fit. This good 
agreement originates from our careful choice of Csw and K in simulations: since we choose these parameters so 
that the region of (M, AM) contains (0, AM^"^) as shown in Fig. ^ csw can be fixed by a short interpolation 
for which the uncertainty due to the choice of the parameterization function for M and AM is not large. 

The alpha's interpolation formula in Ref. (6| gives csw = 2.017 ai (3 — 5.2, which is consistent with our 
results. This confirms that the ALPHA'S formula can be used down to /? = 5.2 as in the UKQCD^IL] and 
JLQCD simulations 12). 

However, as pointed out in Ref. [l9ll. there are large cutoff effects in the PCAC quark mass m and the mass 
dependence of the Sommer scale tq [20j around /9 = 5.2. There is a possibility that the improvement condition 
Eq. H25|l adopted in this and ALPHA'S previous studies leaves unexpectedly large 0{a^) scaling violations in 
physical observables around this value of (3. Therefore, a test of alternative improvement conditions and scaling 
properties of physical observables is an important subject to avoid the large cutoff effects in future lattice 
calculations. 



IV. THREE-FLAVOR QCD 
A. Simulation method 

We determine csw in three-fiavor QCD at nine values of /? in the range /3 = 12.0 -5.2. Numerical simulations 
are carried out on a 8"^ x 16 lattice at four values of csw at each (3, and three or four values of K at each csw- 
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These values are carefully chosen so that the region of (M, AM) contains or is sufficiently close to the point 
(0, AM^°^) which satisfies the improvement condition Eq. (|25|l . These simulation parameters are summarized 
in Table HVI 



In our simulations, we adopt the standard HMC algorithm for two-flavors of dynamical quarks and a polyno- 
mial HMC algorithm developed in Ref. |17| for the remaining one-flavor. We employ the symmetric even-odd 
preconditioning in Refs. (KJIlTl] for the quark matrix D. As in the two-flavor simulations at (3 = 5.2, we calculate 
D~^ using the BiCGStab algorithm with the tolerance parameter ||i?i||/||i?|| < 10^^^. We set the number of 
the Molecular Dynamics steps to 80. This achieves the acceptance rate of about 90% or higher. 

In the PHMC algorithm, we use the Chebyshev polynomial P[D] to approximate D^^. In order to make this 
algorithm exact, the correction factor 

Fcorr = det[W[D]] (33) 

with VF[_D] = P[D]D is taken into account by the noisy Metropolis method [21|. We calculate the square root 
of ly [J], which is required in the Metropolis test, with an accuracy of 10"^"' using the Taylor expansion of 
W[D] [13. The order of the polynomial A^poiy is chosen so that we achieve the acceptance rate of about 90 % 
or higher for the Metropolis test. 

We note that, even with the SF setup, there is a difficulty in simulating massless or negative quark masses 
in three-flavor QCD. In the strong coupling region, eigenvalues of D have large fluctuations and they can take 
values outside the radius of convergence of P[D]. If this happens, the Polynomial approximation P[D] and 
Taylor expansion of M^[_D] break down. For this reason, our simulations in the strong coupling region are 
performed only down to M ~ 0, while negative quark masses M ~ —0.03 are simulated in the weak coupling 
region /3~12. 

We accumulate statistics iVtraj summarized in Table IIVI and measure the correlators fx and /^ {X = A, P) 
at every trajectory. The dependence of the jackknife error of M on the bin size iVbin is investigated in a range 
-^bin = l~-^traj/20. We then adopt A'bin giving the maximum error in the jackknife procedure in the following 
analysis. 

We determine csvv and Kc non-perturbatively also in four-flavor, two-flavor and quenched QCD at /3 = 9.6 to 
study their Nf dependence. The simulation method is similar to that in three-flavor QCD, except that we use 
the standard HMC algorithm in these cases. Simulation parameters are summarized in Tables IVT- IVIII 



B. Non-perturbative csw in three-flavor QCD 

Numerical results of M and AM are summarized in Table IVIIII In Fig. |31 we plot M and M' at several 
values of /? as a function of xq. With our statistics, aM , aM' and hence aAM have an accuracy of 10^"^ at 
all simulation parameters. These accurate data enables us to reduce the statistical error of cgw to the level of 
'^ 5 % even at our coarsest lattice spacing. 

In order to parameterize the K and csw dependence of M and AM, we use the combined fit of Eqs. (|27|l 
and (|28|l . Fit parameters are summarized in Tables IIXI and |X] Figures 0] and El show M and AM and their fit 
curves at several values of /3. We observe that our data are well described by the combined fit. Consequently, as 
shown in Fig.O M dependence of AM is reproduced reasonably well by the fit. We fix csw and Kc satisfying 
the improvement condition Eq. (|25|l using this parameterization. Numerical results for the non-perturbatively 
tuned Csw and Kc are summarized in Table IXll 



As in the analysis in two-flavor QCD at /3 = 5.2, we test Eqs. (|31f) and (|32|) as the alternative method for 
the parametrization. We confirm that the two methods give consistent results both for csw and Kc, and hence 
conclude that the systematic error due to the parameterization method for M and AM is small. 

We fit Csw to a rational function of g^ and obtain the following interpolation formula 



0.194785gg - 0.110781.g;f - 0.0230239gg + 0.137401.gg 
1 - 0.4606855^ 






The interpolation formula for Kc is obtained in a polynomial form 

Kc = 1/8 -ffc(i'<?^ + 0.000964911^0 + 0.00298136ff^-F0.00100995g^-0.00235564g^°, (35) 
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with fc^""^^ =0.00843986. These fits reproduce our data reasonably well with x^/dof of around 1.6. We also note 
that the coefficients oi Olgg) terms in these fits are constrained so that these expressions coincide with their 
one-loop estimates |22,|23| up to O(go)- 

We plot the /3 dependence of csw and Kc in Fig. [7| While csw in three-flavor QCD is well approximated by 
the one-loop estimate in the weak coupling region of Qq <0.4 (/3> 15), it develops a signiflcant deviation toward 
the strong coupling. Similar deviation is also observed in Kc- It is possible that these deviations are partly 
compensated by a better choice of the expansion parameter for lattice perturbation theory [2J|. However, the 
large deviation at the strong coupling (?o ^ 1-0 suggests that one-loop 0(a)-improved Wilson quark action leads 
to a significant 0{a) scaling violation in physical observables at the strong coupling region, where high statistics 
simulations are feasible with currently available computer resources. Therefore, the use of csw hi Eq. H34() is 
essential to remove the 0{a) effects in practical lattice calculations. 



C. Nf dependence 

In Fig. [7| a comparison of our interpolation formula and those by the ALPHA Collaboration in two-flavor and 
quenched QCD suggests that csw monotonously decreases as Nf increases at fixed /3. However, the difference 
between two and three-flavor QCD is not large. This comparison may also suffer from the systematic error due 
to the difference in the analysis method between the two collaborations. In order to study the Nf dependence 
more carefully, we determine csw a-nd Kc for Nf = 4, 2 and at /3 = 9.6 with an analysis method similar to 
that for Nf = 3. We obtai n M and AM summarized in Tables IXlIl - IXIVl Fit to Eqs. ^ and ^ results in 
parameters given in Table IXVI and IXVII Using the improvement condition Eq. H25|l , we obtain csw and Kc 
summarized in Table IXVlII 

Figure IHl shows csw and Kc at (3 = 9.6 as a function of Nf. While Kc has an evident Nf dependence, that for 
Csw is not so clear. The leading Nf dependence of csw is of order g^Nf. By fitting our results to a linear form 
in Nf 

csw = co + ciggiV/, (36) 

we obtain ci = —0.0117(40) which suggests that the two-loop Nf dependence is significant also in csw at the 
relatively weak coupling /3 = 9.6. If two-loop perturbative calculation of csw becomes available in future, it is 
interesting to compare the above estimate of ci to the perturbative estimate. 



D. 0[a/L) uncertainty in csw 

Our non-pcrturbative estimate of csw has 0{a/ L) uncertainties, which affects physical observables at 0{o?) 
when the spatial lattice size L is fixed. However, since we calculate csw with the spatial size in lattice units L/a 
fixed to a constant value 8, the OiajV) dependence of csw induces 0{a) effects in observables. These effects 
can be removed by extrapolating csw to the infinite volume limit, or interpolating to a fixed physical size in 
the whole region of 170. However, we are not able to do that in the present work, since our data are taken at a 
single lattice size at each go- 

We estimate the magnitude of OiajV) uncertainty in csw by using a modified improvement condition. In 
our study as well as in ALPHA'S, the improvement condition Eq. H25|l is adopted instead of AM = in order 
to remove the tree- level OiajV) correction to csw P- ^^ extend this procedure to one- loop level, namely, the 
one- loop correction to AM for the finite lattice volume 8^ x 16 given by |25| 

aAM(i) = -(0.00004839 + 0.000064557V/) (37) 

is incorporated to the improvement condition 

' M =0, 
AM = AM(o)+g2AM(i) ^^^> 

in order to remove the 0{g'^aj V) correction from csw- 

From the parametrization of Eqs. (|27|l - H28|) and the modified improvement condition Eq. H38|) . we obtain 
Csw and Kc summarized in Table IXVIIII A comparison with Table IXII shows that the modified and original 



improvement conditions give consistent results for cgw with each other, and hence that the 0{gQa/ L) correction 
is small in our results. 

In Fig. O we observe that M and Af have a mild xq dependence at Aa < xq < 14a. This may suggest that 
different choices of xq for the improvement conditions Eqs. l|^ and (|S5|l lead to a small difference in csw and 
hence its 0{a/L) ambiguity is not large. 

From these observations, we expect that 0{gQa/L) corrections are not large at L/a — 8, and that csw in the 
infinite volume limit is well approximated by our results. It is an important subject in future studies to confirm 
this point by a direct calculation of csw with varying L. 



V. CONCLUSION 

In this paper, we have performed a non-perturbative 0(a)-improvement of the Wilson quark action in three- 
flavor QCD with the plaquette gauge action. Our high statistics at carefully chosen simulation parameters 
Csw and K enable us to determine non-perturbative csw with an accuracy of ~ 5 % level in the wide range 
oi f]— 12.0-5.2. The main result of this study is the interpolation formula Eq. (|34|) . with which 0{a) scaling 
violation in physical observables can be removed in future simulations at /3 > 5.2. As a byproduct, we also 
obtain the interpolation formula for Kc, which is useful to locate simulation points. 

While it is expected that 0{a/L) uncertainty in csw is not large, this point should be confirmed by a direct 
determination of csw with varying L. A test of scaling properties of physical observables with our estimate of 
Csw is an important next step toward an extensive simulation of three-flavor QCD with the plaquette gauge 
action. 

We note that the range of /3 explored in this study is significantly higher than f3 = 5.0, where an artificial 
phase transition exists. Numerical studies at f3 below our range would possibly suffer from large distortion of 
scaling properties of physical observables. The use of improved gauge actions removes the lattice artifact and 
can push simulations toward coarser lattice spacings. This possibility is explored in a separate publication [Oj. 
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TABLE I: Simulation 


parameters in 


two-flavor QCD at /3=5.2. 






csw 


K 


A^traj 


Csw 


K 


Atraj 


Csw 


K 


At.aj 


1.50 


0.14400 


6400 


1.90 


0.13800 


5600 


2.02 


0.13700 


4800 


1.50 


0.14600 


6200 


1.90 


0.13900 


6387 


2.30 


0.13000 


3200 


1.80 


0.13800 


6400 


2.02 


0.13300 


8000 


2.30 


0.13200 


4000 


1.80 


0.14025 


6400 


2.02 


0.13500 


2176 


3.00 


0.12100 


4000 


1.90 


0.13700 


6400 


2.02 


0.13607 


6400 


3.00 


0.12200 


4000 



TABLE IL Quark mass M and its difference AA/ in two-flavor QCD at /3=5.2. 



csw 



K 



aM 



aAM 



Csw 



K 



aM 



aAM 



Csw 



K 



aM 



aAM 



1.50 0.14400 

1.50 0.14600 

1.80 0.13800 

1.80 0.14025 

1.90 0.13700 



0.159(29) -0.0039(44) 

-0.0012(39) 0.0073(24) 

0.136(10) -0.0073(28) 

-0.0048(29) 0.0024(22) 

0.0629(31) -0.0038(16) 



1.90 0.13800 0.0083(22) -0.0008(16) 

1.90 0.13900 -0.0338(23) 0.0033(14) 

2.02 0.13300 0.1361(51) -0.0003(12) 

2.02 0.13500 0.0541(47) 0.0002(15) 

2.02 0.13607 0.0005(20) -0.0022(14) 



2.02 0.13700 -0.0411(17) 0.0040(27) 

2.30 0.13000 0.0546(10) -0.00108(91) 

2.30 0.13200 -0.0187(11) -0.00296(97) 

3.00 0.12100 0.03767(84) -0.00782(49) 

3.00 0.12200 0.00455(71) -0.00735(82) 





TABLE III: 


Fit parameters for Eqs. 1127^ 


and 


(EHJ in two-flavor QCD at /3 = 


= 5.2. 




x'/dof 


aM 


"m "m 




„(i) 


„{2) 

Cm 




dM 


6.37 


-28.8(7.1) 


8.9(2.5) -0.67(22) 




-6.3(2.2) 


-0.15(17) 




0.83(39) 


xVdof 


aAM 


"am "am 




„(i) 
Cam 


„{2) 

Cam 




rfAM 


1.88 


11.2(3.2) 


-3.9(1.1) 0.336(98) 




3.16(95) 


0.214(71) 




-0.54(17) 



TABLE IV: Simulation parameters in three-flavor QCD. 



11 



/3 = 12.0 



csw = 


= 1.08624 


csw = 


a. 12036 


csw = 


= 1.15448 


csw = 


= 1.18860 


K 


A''traj 


K 


iVtraj 


K 


A''traj 


K 


iVtraj 


0.128958 


3000 


0.128958 


3000 


0.128958 


3000 


0.128958 


3000 


0.129607 


3000 


0.129607 


3000 


0.129607 


3000 


0.129607 


3000 


0.130257 


3000 


0.130257 


3000 


0.130257 


3000 


0.130257 


3000 


0.130906 


3000 


0.130906 


3000 


0.130906 


3000 


0.130906 


3000 



/3 = 9.6 



csw = 


= 1.12990 


csw = 


= 1.16539 


csw = 


= 1.20089 


csw = 


= 1.23638 


K 


A''traj 


K 


iVtraj 


K 


A''traj 


K 


At.aj 


0.13043 


7000 


0.13043 


7000 


0.13043 


7000 


0.13043 


7000 


0.13109 


7000 


0.13109 


7000 


0.13109 


7000 


0.13109 


7000 


0.13175 


7000 


0.13175 


7000 


0.13175 


7000 


0.13175 


7000 


0.13240 


7000 


0.13240 


7000 


0.13230 


1800 


0.13230 


7000 



/3 = 7.4 



csw = 


= 1.2258 


csw = 


= 1.2643 


csw = 


= 1.3028 


csw = 


= 1.3413 


K 


Atraj 


K 


Atraj 


K 


Atraj 


K 


At.aj 


0.13293 


5400 


0.13293 


5400 


0.13293 


5400 


0.13293 


5400 


0.13360 


5400 


0.13360 


5400 


0.13360 


5400 


0.13360 


5400 


0.13427 


5400 


0.13427 


5400 


0.13427 


5400 


0.13393 


400 


0.13494 


5400 


0.13470 


3400 


- 


- 


0.13427 


5400 











/3 = 6.8 








csw = 


= 1.2783 


csw = 


= 1.3184 


csw = 


= 1.3586 


csw = 


= 1.3987 


K 


Atraj 


K 


Atraj 


K 


Atraj 


K 


At.aj 


0.13391 


4200 


0.13391 


4200 


0.13391 


4200 


0.13391 


4200 


0.13459 


4200 


0.13459 


4200 


0.13459 


4200 


0.13459 


4200 


0.13526 


4200 


0.13526 


4200 


0.13526 


4200 


0.13492 


3700 


0.13594 


4200 


0.13560 


4200 


0.13540 


3700 


0.13500 


3150 











/3 = 6.3 








csw = 


= 1.3117 


csw = 


= 1.3675 


csw = 


= 1.4233 


csw = 


= 1.4791 


K 


At.aj 


K 


At.aj 


K 


At.aj 


K 


At.aj 


0.13501 


4600 


0.13501 


4600 


0.13446 


4600 


0.13446 


4600 


0.13555 


4600 


0.13555 


4600 


0.13501 


4600 


0.13501 


4600 


0.13609 


4600 


0.13609 


4600 


0.13555 


4600 


0.13530 


1100 


0.13664 


4600 


0.13664 


4600 


0.13609 


4600 


0.13555 


4600 











/3 = 6.0 








csw = 


= 1.3237 


csw = 


= 1.3801 


csw = 


= 1.4364 


csw = 


= 1.4927 


K 


At.aj 


K 


At.aj 


K 


At.aj 


K 


At.aj 


0.13552 


4200 


0.13552 


4200 


0.13498 


4200 


0.13498 


4200 


0.13607 


4200 


0.13607 


4200 


0.13552 


4200 


0.13552 


4200 


0.13661 


4200 


0.13661 


4200 


0.13607 


4200 


0.13607 


4200 


0.13716 


4200 


0.13716 


4200 


0.13661 


4200 


0.13634 


4200 











/3 = 5.7 








csw = 


= 1.5431 


csw = 


= 1.6088 


csw = 


= 1.6745 


csw = 


= 1.7401 


K 


Atraj 


K 


Atraj 


K 


Atraj 


K 


Atraj 


0.13575 


5000 


0.13541 


5000 


0.13507 


5000 


0.13400 


4500 


0.13644 


5000 


0.13610 


5000 


0.13540 


4400 


0.13450 


4100 


0.13712 


4700 


0.13644 


3700 


0.13575 


5000 


0.13473 


5000 


0.13750 


4500 


0.13678 


5000 


0.13585 


2700 


0.13500 


4400 
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TABLE IV; (Continued) 



/3 = 5.4 



csw- 


= 1.6 


csw = 


= 1.7 


csw = 


= 1.8 


csw = 


= 1.9 


K 


iVtraj 


K 


iVt.aj 


K 


iVt.aj 


K 


iVtraj 


0.13750 


9000 


0.13600 


9000 


0.13480 


9000 


0.13330 


9000 


0.13810 


9000 


0.13630 


7900 


0.13520 


9000 


0.13370 


9000 


0.13825 


7100 


0.13660 


9000 


0.13560 


9000 


0.13410 


9000 


0.13840 


8200 


0.13720 


9000 


0.13580 


7700 


0.13450 


9000 



13 = 5.2 



csw = 


= 1.70 


csw = 


= 1.85 


csw - 


= 2.00 


csw = 


= 2.15 


K 


iVtraj 


K 


A^traj 


K 


Atraj 


K 


At.aj 


0.1375 


7800 


0.1357 


7800 


0.1336 


7800 


0.1316 


7800 


0.1379 


7800 


0.1361 


7800 


0.1340 


7800 


0.1320 


7800 


0.1383 


7800 


0.1365 


7800 


0.1344 


7800 


0.1322 


5500 


0.1387 


6300 


0.1366 


7400 


0.1345 


7400 


0.1324 


7800 



TABLE V: Simulation parameters in four-flavor QCD at /3 = 9.6. 



csw = 


= 1.12990 


csw = 


= 1.16539 


csw = 


= 1.20089 


csw = 


= 1.23638 


K 


Atraj 


K 


Atraj 


K 


At.aj 


K 


Atraj 


0.13043 


5000 


0.13043 


5000 


0.13043 


5000 


0.13043 


5000 


0.13109 


5000 


0.13109 


5000 


0.13109 


5000 


0.13109 


5000 


0.13175 


5000 


0.13175 


5000 


0.13175 


5000 


0.13175 


5000 


0.13240 


5000 


0.13240 


5000 


0.13240 


5000 


0.13240 


5000 



TABLE VI: Simulation parameters in two-flavor QCD at /3 = 9.6. 



csw = 


= 1.12990 


csw = 


= 1.16539 


csw = 


= 1.20089 


csw = 


= 1.23638 


K 


At.aj 


K 


At.aj 


K 


At.aj 


K 


At.aj 


0.13043 


5000 


0.13043 


5000 


0.13043 


5000 


0.13043 


5000 


0.13109 


5000 


0.13109 


5000 


0.13109 


5000 


0.13109 


5000 


0.13175 


5000 


0.13175 


5000 


0.13175 


5000 


0.13175 


5000 


0.13240 


5000 


0.13240 


5000 


0.13240 


5000 


0.13240 


5000 



TABLE VII: Simulation parameters in quenched QCD at /3 = 9.6. 



csw = 


= 1.12990 


csw = 


= 1.16539 


csw = 


= 1.20089 


csw = 


= 1.23638 


K 


Atraj 


K 


Atraj 


K 


Atraj 


K 


Atraj 


0.13043 


8000 


0.13043 


8000 


0.13043 


8000 


0.13043 


8000 


0.13109 


8000 


0.13109 


8000 


0.13109 


8000 


0.13109 


8000 


0.13175 


8000 


0.13175 


8000 


0.13175 


8000 


0.13175 


8000 


0.13240 


8000 


0.13240 


8000 


0.13240 


8000 


0.13240 


8000 



TABLE VIII: Quark mass M and its difference AM in three-flavor QCD. 
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/3 = 12.0 



csw = 1.08624 


csw = 1.12036 


K 


aM 


aAM 


K 


aM 


aAM 


0.128958 


0.03502(11) 


0.00107(13) 


0.128958 


0.03066(8) 


0.00067(10) 


0.129607 


0.01498(13) 


0.00099(12) 


0.129607 


0.01028(10) 


0.00072(16) 


0.130257 


-0.00544(10) 


0.00100(14) 


0.130257 


-0.01011(12) 


0.00040(12) 


0.130906 


-0.02591(11) 


0.00116(12) 


0.130906 


-0.03088(11) 


0.00052(12) 


csw = 1.15448 


csw = 1.18860 


K 


aM 


aAM 


K 


aM 


aAM 


0.128958 


0.02581(9) 


0.00027(12) 


0.128958 


0.02125(9) 


-0.00029(11) 


0.129607 


0.00561(11) 


0.00015(11) 


0.129607 


0.00076(10) 


-0.00027(11) 


0.130257 


-0.01512(14) 


0.00008(13) 


0.130257 


-0.01965(11) 


-0.00044(13) 


0.130906 


-0.03555(17) 


0.00025(13) 


0.130906 


-0.04041(12) 


-0.00043(14) 



/3 = 9.6 



csw = 1.12990 


csw = 1.16539 


K 


aM 


aAAf 


K 


aM 


aAM 


0.13043 


0.03818(12) 


0.00118(10) 


0.13043 


0.03226(9) 


0.00067(10) 


0.13109 


0.01805(9) 


0.00100(12) 


0.13109 


0.01187(10) 


0.00062(12) 


0.13175 


-0.00275(10) 


0.00104(14) 


0.13175 


-0.00871(12) 


0.00056(13) 


0.13240 


-0.02312(11) 


0.00104(13) 


0.13240 


-0.02948(11) 


0.00063(16) 


csw = 1.20089 


csw = 1.23638 


K 


aM 


aAM 


K 


aM 


aAM 


0.13043 


0.02596(11) 


0.00037(11) 


0.13043 


0.01995(9) 


-0.00037(14) 


0.13109 


0.00559(12) 


0.00014(13) 


0.13109 


-0.00091(10) 


-0.00012(9) 


0.13175 


-0.01522(11) 


0.00005(12) 


0.13175 


-0.02145 (11) 


-0.00044(13) 


0.13230 


-0.03256(19) 


0.00033(24) 


0.13230 


-0.03934 (13) 


-0.00054(14) 



/3 = 7.4 



csw = 1.2258 


csw = 1.2643 


K 


aM 


aAM 


K 


aM 


aAM 


0.13293 


0.04147(22) 


0.00135(19) 


0.13293 


0.03265(15) 


0.00067(19) 


0.13360 


0.02134(21) 


0.00158(25) 


0.13360 


0.01209(23) 


0.00065(25) 


0.13427 


0.00079(21) 


0.00099(26) 


0.13427 


-0.00942(21) 


0.00102(23) 


0.13494 


-0.02065(24) 


0.00077(28) 


0.13470 


-0.02292(28) 


0.00089(34) 


csw = 1.3028 


csw = 1.3413 


K 


aM 


aAM 


K 


aM 


aAM 


0.13293 


0.02296(25) 


0.00035(24) 


0.13293 


0.01309(22) 


0.00004(23) 


0.13360 


0.00214(17) 


0.00018(19) 


0.13360 


-0.00756(21) 


0.00029(28) 


0.13427 


-0.01888(23) 


0.00052(22) 


0.13393 


-0.01774(64) 


-0.00018(68) 


- 


- 


- 


0.13427 


-0.02867(26) 


-0.00048(36) 



/3 = 6.8 



csw = 1.2783 


csw = 1.3184 


K 


aM 


aAM 


K 


aM 


aAM 


0.13391 
0.13459 
0.13526 
0.13594 


0.04462(36) 

0.02457(29) 

0.00338(33) 

-0.01824(30) 


0.00066(28) 
0.00079(32) 
0.00098(30) 
0.00163(37) 


0.13391 
0.13459 
0.13526 
0.13560 


0.03342(33) 

0.01358(28) 

-0.00823(30) 

-0.01839(36) 


0.00092(28) 
0.00051(27) 
0.00065(40) 
0.00054(39) 


csw = 1.3586 


csw = 1.3987 


K 


aM 


aAM 


K 


aM 


aAM 


0.13391 
0.13459 
0.13526 
0.13540 


0.02240(31) 

0.00175(31) 

-0.01924(27) 

-0.02474(30) 


0.00049(26) 
0.00031(23) 
0.00033(34) 
0.00058(46) 


0.13391 
0.13459 
0.13492 
0.13500 


0.01208(22) 
-0.00953(30) 
-0.01956(58) 
-0.02271(43) 


-0.00021(35) 

-0.00030(27) 

0.00025(45) 

-0.00014(34) 



TABLE VIII; (Continued) 
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/3 = 6.3 



csw = 1.3117 


csw = 1.3675 


K 


aM 


aAM 


K 


aM 


aAM 


0.13501 


0.05392(55) 


0.00085(50) 


0.13501 


0.03677(44) 


0.00097(29) 


0.13555 


0.03766(40) 


0.00107(32) 


0.13555 


0.02065(30) 


0.00063(33) 


0.13609 


0.02152(36) 


0.00080(38) 


0.13609 


0.00353(36) 


0.00058(32) 


0.13664 


0.00403(41) 


0.00175(34) 


0.13664 


-0.01345(53) 


0.00168(41) 


csw = 1.4233 


csw = 1.4791 


K 


aM 


aAM 


K 


aM 


aAM 


0.13446 


0.03635(37) 


0.00070(37) 


0.13446 


0.01834(37) 


-0.00034(33) 


0.13501 


0.01915(34) 


0.00028(37) 


0.13501 


0.00198(38) 


0.00006(30) 


0.13555 


0.00284(43) 


0.00092(33) 


0.13530 


-0.00807(57) 


0.00075(47) 


0.13609 


-0.01497(43) 


0.00050(37) 


0.13555 


-0.01599(41) 


-0.00039(51) 









13 = 6.0 






csw = 1.3237 


csw = 1.3801 


K 


aM 


aAM 


K 


aM 


aAM 


0.13552 
0.13607 
0.13661 
0.13716 


0.07608(51) 
0.06056(57) 
0.04411(52) 
0.02755(53) 


0.00154(39) 
0.00175(43) 
0.00120(48) 
0.00241(53) 


0.13552 
0.13607 
0.13661 
0.13716 


0.05826(50) 
0.04182(39) 
0.02400(46) 
0.00843(50) 


0.00156(37) 
0.00124(44) 
0.00077(40) 
0.00259(53) 


csw = 1.4364 


csw = 1.4927 


K 


aM 


aAM 


K 


aM 


aAM 


0.13498 
0.13552 
0.13607 
0.13661 


0.05492(48) 
0.03893(49) 
0.02221(48) 
0.00562(49) 


0.00054(40) 
0.00035(54) 
0.00107(44) 
0.00031(50) 


0.13498 
0.13552 
0.13607 
0.13634 


0.03635(37) 

0.01878(41) 

0.00235(45) 

-0.00692(49) 


-0.00060(60) 
0.00030(45) 

-0.00005(48) 
0.00076(39) 



13 = 5.7 



csw = 1.5431 


csw = 1.6088 


K 


aM 


aAM 


K 


aM 


aAM 


0.13575 


0.04354(59) 


0.00079(42) 


0.13541 


0.02748(50) 


-0.00068(47) 


0.13644 


0.02091(55) 


-0.00039(51) 


0.13610 


0.00695(130) 


0.00053(54) 


0.13712 


-0.00181(57) 


0.00072(53) 


0.13644 


-0.00587(70) 


-0.00087(64) 


0.13750 


-0.01386(70) 


0.00007(72) 


0.13678 


-0.01702(60) 


-0.00089(65) 


csw = 1.6745 


csw = 1.7401 


K 


aM 


aAM 


K 


aM 


aAM 


0.13507 


0.01331(69) 


-0.00081(51) 


0.13400 


0.02229(42) 


-0.00140(52) 


0.13540 


0.00278(64) 


-0.00022(61) 


0.13450 


0.00521(59) 


-0.00104(39) 


0.13575 


-0.00947(54) 


-0.00078(43) 


0.13473 


-0.00173(44) 


-0.00082(42) 


0.13585 


-0.01124(62) 


-0.00115(51) 


0.13500 


-0.01153(59) 


-0.00172(69) 



13 = 5.4 



csw = 1.6 


csw = 1.7 


K 


aM 


aAM 


K 


aM 


aAM 


0.13750 
0.13810 
0.13825 
0.13840 


0.03310(70) 
0.01245(75) 
0.00684(78) 
0.00382(124) 


0.00015(43) 
0.00083(52) 
0.00053(62) 
0.00100(49) 


0.13600 
0.13630 
0.13660 
0.13720 


0.03592(56) 

0.02622(54) 

0.01653(48) 

-0.00589(52) 


-0.00031(40) 
0.00032(52) 
0.00030(46) 
0.00106(57) 


csw = 1.8 


csw = 1.9 


K 


aM 


aAM 


K 


aM 


aAM 


0.13480 
0.13520 
0.13560 
0.13580 


0.03073(51) 

0.01650(50) 

0.00373(52) 

-0.00426(70) 


-0.00129(56) 

-0.00048(49) 

0.00010(47) 

-0.00062(53) 


0.13330 
0.13370 
0.13410 
0.13450 


0.03600(40) 

0.02250(34) 

0.00871(47) 

-0.00508(51) 


-0.00149(40) 
-0.00167(51) 
-0.00110(40) 
-0.00194(51) 



TABLE VIII: {Continued) 



(3 = 5.2 
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csw = 1.70 


csw = 1.85 


K 


aM 


aAM 


K 


aM 


aAM 


0.1375 
0.1379 
0.1383 
0.1387 


0.05191(131) 
0.03776(89) 
0.02311(82) 
0.00807(123) 


0.00102(87) 
0.00062(67) 
0.00196(54) 
0.00040(110) 


0.1357 
0.1361 
0.1365 
0.1366 


0.03406(87) 
0.01844(73) 
0.00287(78) 
0.00037(84) 


-0.00004(54) 
0.00045(46) 
0.00019(81) 

-0.00038(56) 


csw = 2.00 


csw = 2.15 


K 


aM 


aAM 


K 


aM 


aAM 


0.1336 
0.1340 
0.1344 
0.1345 


0.02982(64) 

0.01467(74) 

-0.00003(68) 

-0.00405(69) 


-0.00161(43) 
-0.00089(48) 
-0.00219(48) 
-0.00133(48) 


0.1316 
0.1320 
0.1322 
0.1324 


0.02568(70) 

0.01350(49) 

0.00364(74) 

-0.00202(51) 


-0.00314(65) 
-0.00281(46) 
-0.00348(77) 
-0.00284(60) 



TABLE IX: Fit parameters for Eq. ^ in throo-flavor QCD. 



xVdof 



Um 



^(1) 



k(2) 



.(1) 



„(2) 



12.0 


0.96 


-10.9( 1.1) 


2.40(28) 


-0.125(18) 


-0.47(13) 


-0.012(23) 


0.046(16) 


9.6 


1.41 


-13.1(1.1) 


3.02(30) 


-0.168(20) 


-0.64(14) 


-0.047(22) 


0.076(17) 


7.4 


1.17 


-15.9(2.9) 


3.90(81) 


-0.233(55) 


-0.90(34) 


-0.057(39) 


0.106(46) 


6.8 


1.57 


-24.9(4.3) 


6.6(1.2) 


-0.430(85) 


-2.30(47) 


0.016(53) 


0.266(68) 


6.3 


0.66 


-19.0(5.8) 


5.0(1.7) 


-0.33(12) 


-1.99(66) 


-0.094(40) 


0.262(97) 


6.0 


1.19 


-12.2(6.6) 


3.2(1.9) 


-0.20(13) 


-1.75(64) 


-0.111(40) 


0.234(92) 


5.7 


1.16 


-7.6(9.5) 


1.6(3.0) 


-0.07(23) 


-0.4(1.8) 


-0.007(99) 


0.005(280) 


5.4 


0.92 


-27(11) 


8.0(3.7) 


-0.58(30) 


-5.0(2.7) 


-0.19(17) 


0.70(45) 


5.2 


1.42 


-38(22) 


11.9(7.4) 


-0.91(62) 


-8.0(5.5) 


-0.32(35) 


1.17(93) 



TABLE X: Fit parameters for Eq. ^ in three-flavor QCD. 



xVdof 



aAM 



hW 



!,(2) 



„(1) 



„(2) 



dAJ\ 



12.0 


0.59 


1.6( 1.2) 


-0.41(32) 


0.025(21) 


-0.12(16) 


-0.008(26) 


0.016(19) 


9.6 


0.91 


0.6(1.3) 


-0.14(35) 


0.008(23) 


-0.14(16) 


-0.007(25) 


0.019(19) 


7.4 


1.13 


-3.1(3.4) 


0.82(93) 


-0.053(64) 


0.09(38) 


0.025(44) 


-0.022(51) 


6.8 


0.40 


6.2(4.4) 


-1.7(1.2) 


0.109(84) 


-0.06(50) 


-0.051(49) 


0.025(66) 


6.3 


1.14 


7.1(5.3) 


-2.0(1.5) 


0.14(11) 


0.49(61) 


-0.012(38) 


-0.062(91) 


6.0 


1.28 


8.7(6.3) 


-2.5(1.8) 


0.18(13) 


0.78(60) 


0.003(37) 


-0.108(84) 


5.7 


1.06 


-4.8(8.7) 


1.5(2.7) 


-0.11(21) 


-0.7(1.5) 


0.008(80) 


0.09(23) 


5.4 


0.53 


-6.9(10.8) 


2.4(3.6) 


-0.21(29) 


-2.1(2.6) 


-0.18(16) 


0.38(43) 


5.2 


0.67 


-26(16) 


8.7(5.4) 


-0.73(45) 


-6.6(4.0) 


-0.42(25) 


1.10(67) 



TABLE XI: Non-perturbative estimate of csw and Kc in three-flavor QCD obtained from tree-level improvement 
condition Eq. 11251 . 



P 



Csw 



K, 



12.0 

9.6 

7.4 

6.8 

6.3 

6.0 

5.7 

5.4 

5.2 



1.1415(48) 
1.1916(50) 
1.316(11) 
1.358(13) 

1.447(15) 
1.494(14) 
1.544(32) 
1.740(30) 
1.764(103) 



0.129841(21) 

0.131321(28) 

0.133567(87) 

0.13466(12) 

0.13539(16) 

0.13612(16) 

0.13706(39) 

0.13650(41) 

0.13789(159) 
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TABLE XII: Quark mass M and its difference AA/ in four-flavor QCD at /3 = 9.6. 



csw = 1.12990 


csw = 1.16539 


K 


aM 


aAM 


K 


aM 


aAM 


0.13043 
0.13109 
0.13175 
0.13240 


0.03530(11) 

0.01527(12) 

-0.00522(12) 

-0.02611(13) 


0.00111(14) 
0.00118(15) 
0.00110(12) 
0.00123(12) 


0.13043 
0.13109 
0.13175 
0.13240 


0.02955(11) 

0.00892(12) 

-0.01158(11) 

-0.03256(14) 


0.00050(11) 
0.00068(13) 
0.00069(12) 
0.00061(13) 


csw = 1.20089 


csw = 1.23638 


K 


aM 


aAAf 


K 


aM 


aAA/ 


0.13043 
0.13109 
0.13175 
0.13240 


0.02328(12) 

0.00280(14) 

-0.01825(10) 

-0.03903(15) 


0.00018(11) 
0.00028(13) 
0.00028(15) 
0.00004(15) 


0.13043 
0.13109 
0.13175 
0.13240 


0.01723(12) 
-0.00364(12) 
-0.02462(17) 
-0.04592(14) 


-0.00032(14) 
-0.00021(11) 
-0.00035(15) 
-0.00023(13) 



TABLE XIII: Quark mass M and its difference AM in two-flavor QCD at /3 = 9.6. 



csw = 1.12990 


csw = 1-16539 


K 


aM 


aAA'f 


K 


aM 


aAM 


0.13043 
0.13109 
0.13175 
0.13240 


0.04089(10) 

0.02085(10) 

0.00047(12) 

-0.01998(13) 


0.00127(13) 
0.00119(13) 
0.00102(16) 
0.00104(13) 


0.13043 
0.13109 
0.13175 
0.13240 


0.03508(10) 

0.01482(12) 

-0.00585(13) 

-0.02646(15) 


0.00085(12) 
0.00053(16) 
0.00070(14) 
0.00068(15) 


csw = 1.20089 


csw = 1.23638 


K 


aM 


aAAf 


K 


aM 


aAA/ 


0.13043 
0.13109 
0.13175 
0.13240 


0.02884(12) 

0.00855(11) 

-0.01225(13) 

-0.03270(12) 


0.00013(13) 
0.00045(14) 
0.00013(19) 
0.00032(22) 


0.13043 
0.13109 
0.13175 
0.13240 


0.02259(10) 

0.00218(10) 

-0.01857(13) 

-0.03927(11) 


-0.00021(13) 
-0.00003(12) 
-0.00000(13) 
-0.00024(16) 



TABLE XIV: Quark mass M and its difference AM in quenched QCD at /3 = 9.6. 



csw = 1.12990 


csw = 1.16539 


K 


aM 


aAM 


K 


aM 


aAA/ 


0.13043 
0.13109 
0.13175 
0.13240 


0.04646(7) 

0.02679(8) 

0.00656(7) 

-0.01324(8) 


0.00141(10) 
0.00131(11) 
0.00131(14) 
0.00132(10) 


0.13043 
0.13109 
0.13175 
0.13240 


0.04059(8) 

0.02064(8) 

0.00041(7) 

-0.01978(8) 


0.00095(11) 
0.00080(11) 
0.00084(10) 
0.00074(11) 


csw = 1.20089 


csw = 1.23638 


K 


aM 


aAAf 


K 


aAf 


aAM 


0.13043 
0.13109 
0.13175 
0.13240 


0.03458(9) 

0.01451(7) 

-0.00580(8) 

-0.02607(10) 


0.00036(11) 
0.00034(12) 
0.00050(14) 
0.00053(16) 


0.13043 
0.13109 
0.13175 
0.13240 


0.02840(8) 

0.00818(7) 

-0.01231(8) 

-0.03247(10) 


-0.00015(10) 

0.00006(10) 

-0.00005(13) 

-0.00013(21) 









TABLE XV: 


Fit 


parameters 


in Eq. (HTJ for Nf 


= 4, 2 and 0. 






Nf 


P 


xVdof 


aM 




"m 


"m 


„(i) 


Cm 


dM 


4 


9.6 


0.96 


-14.7(1.3) 




3.52(33) 


-0.206(22) 


-1.14(15) 


-0.041(25) 


0.140(19) 


2 


9.6 


0.48 


-13.5(1.2) 




3.15(31) 


-0.177(20) 


-0.64(14) 


-0.052(23) 


0.077(17) 





9.6 


1.03 


-11.71(80) 




2.70(21) 


-0.150(14) 


-0.80(10) 


-0.040(16) 


0.094(12) 
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TABLE XVI: Fit parameters in Eq. ^ for iV/ =4, 2 and 0. 



Nf 


P 


xVdof 


a AM 


"am 


"am 


'-AM 


'-AM 


dAM 


4 


9.6 


0.42 


-1.2(1.3) 


0.33(34) 


-0.023(22) 


-0.14(16) 


0.015(26) 


0.012(19) 


2 


9.6 


0.80 


-1.0(1.4) 


0.24(38) 


-0.014(25) 


0.15(17) 


0.013(28) 


-0.026(21) 





9.6 


0.53 


-0.3(1.2) 


0.06(31) 


-0.002(20) 


0.11(14) 


0.012(23) 


-0.020(18) 



TABLE XVII: Non-perturbative estimate of csw and Kc for Nf = 4, 2 and obtained from tree-level improvement 
condition Eq. 1251 . 



Nf 


P 


Csw 


K, 


4 


9.6 


1.1954(48) 


0.131210(27) 


2 


9.6 


1.2028(63) 


0.131353(36) 





9.6 


1.2112(44) 


0.131502(25) 



TABLE XVIII: Non-perturbative estimate of csw and Kc in three-flavor QCD obtained from one-loop level improvement 
condition Eq. (PHJ- 



/3 



csw 



K, 



12.0 

9.6 

7.4 

6.8 

6.3 

6.0 

5.7 

5.4 

5.2 



1.1507(46) 

1.2039(46) 

1.336(14) 

1.378(12) 

1.470(16) 

1.508(17) 

1.569(27) 

1.770(26) 

1.811(67) 



0.129801(20) 

0.131252(26) 

0.13340(12) 

0.13447(11) 

0.13515(17) 

0.13596(20) 

0.13676(32) 

0.13609(35) 

0.1372(10) 
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FIG. 1: Plots of M (top figure) and AM (bottom figure) in two-flavor QCD at /3 = 5.2 as a function of 1/K. The 
bottom figure shows AM as a function of M. 
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N=2, 8^16, B=5.2 
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FIG. 2: Mass difference AM at M = as a function of csw. Filled and open symbols are obtained from the fit Eq. I{'M^ 
and combined fit Eqs. 1271 and 12811 . Results from the combined fit are slightly shifted in the horizontal direction for a 
better visibility. Dotted and dashed lines show linear and quadratic fit Eq. 13211 . Dot-dashed line shows the tree- level 
value AM™ =0.000277. 
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FIG. 3: Quark masses M and M' in three-flavor QCD as a function of xo. Top, middle and bottom flgures show data 
at /3=12.0, 9.6 and 5.2, respectively. 
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FIG. 4: Quark mass M as a function of 1/K. Top, middle and bottom figures show data at /3 = 12.0, 9.6 and 5.2, 
respectively. Dotted lines are fit curve of Eq. 12711 . 
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FIG. 5: Mass difference AM as a function of 1/K. Top, middle and bottom figures show data at /3=12.0, 9.6 and 5.2, 
respectively. Dotted lines are fit curve of Eq. H28II . 
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FIG. 6: Mass difference AM as a function of M in three-flavor QCD. Top, middle and bottom figures show data at 
(3=12.0, 9.6 and 5.2, respectively. Dotted lines are fit curve reproduced from Eqs. 127^ and 12811 . 
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FIG. 7: Non-perturbatively determined csw (top figure) and Kc (bottom figure) as a function of go- For Kc in quencfied 
QCD, autlior's interpolation of tlie ALPHA'S results is plotted. 
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FIG. 8: Plots of csw (top figure), Kc (bottom figure) at /3 = 9.6 as a function of Nf. Dashed line shows linear fit to 
data. 



